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ABSTRACT 

We present a new approach to measure the mass function of dark matter halos and to discriminate 
models with differing values of f2 through weak gravitational lensing. We measure the distribution of 
peaks from simulated lensing surveys and show that the lensing signal due to dark matter halos can 
be detected for a wide range of peak heights. Even when the signal-to-noise is well below the limit 
for detection of individual halos, projected halo statistics can be constrained for halo masses spanning 
galactic to cluster halos. The use of peak statistics relies on an analytical model of the noise due to the 
intrinsic ellipticities of source galaxies. The noise model has been shown to accurately describe simulated 
data for a variety of input ellipticity distributions. We show that the measured peak distribution has 
distinct signatures of gravitational lensing, and its non-Gaussian shape can be used to distinguish models 
with different values of 17. The use of peak statistics is complementary to the measurement of field 
statistics, such as the ellipticity correlation function, and possibly not susceptible to the same systematic 
errors. 
Subject headings: cosmology: theory — cosmology: gravitational lensing — methods: numerical 



1. INTRODUCTION 

In the coming years, gravitational lensing is likely to 
become an effective tool for mapping large-scale structure 
in the universe. Over the past decade several measure- 
ments of weak lensing by galaxy clusters have been made. 
Mass reconstruction techniques are now being applied to 
wide field lensing surveys in blank fields that will probe 
the dark matter distribution over angular scales of order 
1' — 1°. Wide field lensing observations have already de- 
tected filaments and dark halos that were not visible by 
their light distribution (Kaiser at al 1998; Erben et al 1999; 
Tyson et al. 1999). 

Statistical properties of the clustering of dark matter 
can be probed from lensing data by computing shear cor- 
relations over blank fields with area of order 10 square de- 
grees (Blandford et al 1991; Miralda-Escude 1991; Kaiser 
1992; Bernardeau et al 1997; Jain & Seljak 1997; Kaiser 
1998; Stebbins 1996; Schneider et al. 1998). An alterna- 
tive approach is to focus on the statistics of dark matter 
halos, identified through their lensing strength, using mea- 
sures such as the aperture mass (Schneider 1996; Kruse & 
Schneider 1999a; Kruse & Schneider 1999a; Reblinsky et 
al. 1999). The halo statistics approach has been shown by 
the above authors to be a useful probe of the mass func- 
tion for massive, cluster sized halos; the main practical 
limitation is that only ~ 10 halos per square degree are 
expected to be detected with adequate signal-to-noise. 

This paper advocates a new approach to the measure- 
ment of the statistics of dark matter halos through lensing. 
By modeling the distribution of peaks in lensing data in- 
duced by the noise due to the intrinsic ellipticities of source 
galaxies, we show that it is possible to statistically detect 
the signal due to dark matter halos, even for mass scales 
below the signal-to-noise limit for the detection of indi- 
vidual halos. Section 2 describes the construction of peak 
statistics from simulated data and from pure noise. Re- 



sults for the peak statistics for a set of cosmological models 
are shown in Section 3. We discuss the prospects for mea- 
suring the halo mass function and discriminating models 
from realistic data in Section 4. 

2. PEAK STATISTICS IN SIMULATED DATA 

We use shear and convergence fields from ray tracing 
simulations through the dark matter distribution of N- 
body simulations (Jain, Seljak & White 1999). The fields 
we use are about 3 degrees on a side, sampled with a 
grid spacing of 0.1', with source galaxies taken to be at 
z — 1. We use two cosmological models, an Einstein-de 
Sitter model and an open model with r^mattcr = 0.3. The 
power spectrum corresponds to a cold dark matter shape 
parameter F = 0.21 model. Further details of the models 
and the simulations are given in Jain et al (1999). 

A simulated noisy map of the convergence, k,{9), is built 
by first smoothing the k field over scale 9q with a Gaussian 
window W{e) = e-x:p{-\e\^ /el;)/Ti0%. The noise due to the 
randomly oriented intrinsic ellipticities of source galaxies 
is modeled as a Gaussian random field with variance, 
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where a^ is the rms amplitude of the intrinsic ellipticity 
distribution and ng is the number density of source galax- 
ies. This Gaussian noise is added to the smoothed k field; 
the accuracy of this noise model is discussed below. From 
the smoothed noisy data, peaks are found by identifying 
pixels that have a higher/lower value of k than all neigh- 
boring pixels. This corresponds to the condition that the 
gradient of the field vanishes and thus includes peaks as 
well as troughs. The height of the peak r/ = K/tTnoise is its 
value in units of the noise rms in the smoothed field. 

Our choice of the noise model for the convergence field 
relies on previous work. Van Waerbeke et al (1999) have 
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Fig. 1. — Histogram of the peak distribution in a noisy convergence field for the open model. The solid line with error bars is the peak 
distribution measured from the k map with Gaussian noise added as discussed in the text. The dashed line is the peak distribution measured 
from the k map reconstructed from the ellipticities. The near coincidence of the two curves demonstrates the accuracy of the reconstruction 
scheme and of the noise model. The error bars are computed from 7 realizations of the signal k map. The dot-dashed line shows the peak 
distribution due to the noise alone. Almost overlapping with it is the analytical model for noise peaks discussed in the text. 




Fig. 2. — Comparison of the reconstructed and input maps of k. The left panel shows the input k maps for a field 3 degrees on a side, 
and the right panel shows the same field reconstructed from cllipticity data that includes the noise due to the intrinsic ellipticities of galaxies. 
While the small scale peaks in the right panel are dominated by the noise, the subsequent figures show how their distribution contains imprints 
of the signal. 



shown that the convergence field can be accurately re- 
constructed from observed ellipticity data in the absence 
of systematic errors. Both the reconstruction schemes of 
Kaiser & Squires (1993) and the Maximum-Likelihood al- 
gorithm of Bartelmann et al (1996) recover the conver- 
gence field with adequate accuracy for fields of order a de- 
gree on a side. Van Waerbeke (1999) further showed that 
the noise properties of peaks can be analytically described 
using Gaussian statistics (Bardeen et al 1986; Bond & Ef- 
stathiou 1987), and the weak lensing approximation. 

Figure |l| shows a test of the analytical model of Van 
Waerbeke (1999) for peaks due to noise, and checks the 
accuracy of the peak distribution in the reconstructed n. 
The dot-dashed curve shows the histogram of peaks from 
this analytical noise model; the double peaked shape is due 
to peaks with positive curvature and troughs with negative 
curvature. Almost overlapping with the analytical curve 
is the measured histograms of peaks in a field with pure 
noise. The distribution of peak heights measured in maps 
of K reconstructed from noisy ellipticity data (dashed line) 



is compared with the distribution in the k maps of the si, 
nal plus Gaussian noise with variance given by equation 
(solid line). The close agreement of the two curves demon- 
strates the accuracy of the k reconstruction scheme. 

The success of the reconstruction gives us confidence in 
working with the convergence data and the noise model 
directly, avoiding the slow and expensive reconstruction 
process. We have also verified that using the mass aper- 
ture statistic, which is constructed directly from ellipticity 
data, leads to very similar peak distributions. We have 
also found that a variety of distributions of the intrinsic 
ellipticity (including non-Gaussian ditributions) produce 
the same Gaussian statistics for peaks in the n maps; these 
results will be presented elsewhere. Figure g shows the ac- 
tual maps of the convergence, k, used to measure the peak 
distributions. The small amplitude peaks in the signal 
map are swamped by the noise, so there is little hope of 
recovering them individually from data. However, we show 
below that their distribution is sufficiently modulated by 
the signal to distinguish cosmological models. 
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Fig. 3. — Probability distribution function of the peaks in noise-free fields (left panels) and noisy fields (right panels) 3 degrees on a side, 
with smoothing scale 9g = 0.5' (upper panels) and 1' (lower panels). The solid lines show the open model and the dashed lines the EdS 
model. The error bars are obtained from averaging over seven realizations. In the right panels the dot-dashed lines show the pdf of peaks in 
fields with pure noise. 



3. SENSITIVITY OF PEAK STATISTICS TO Q, 

Figure y shows the probabiHty distribution function 
(pdf) of peaks in the convergence field as a function of peak 
height for noise-free (left panels) and noisy (right panels) 
fields for two different smoothing scales. The pdf in noise- 
free fields has the qualitative characteristics due to non- 
linear gravitational clustering: at negative peak heights 
(underdense regions) it has a cutoff related to the mini- 
mum K resulting from empty beams, and it has a tail at 
positive v due to collapsed halos. The cosmological mod- 
els have different pdf 's, just as they do for the pdf of k in 
the field, shown in figure Q. For the noisy fields, the num- 
ber density of source galaxies is 30 per square arcminute, 
and their rms intrinsic ellipticity is CTg = 0.2. Thus the 
peak height v = 1 corresponds to an averaged value of 
K = 0.02 (O.Of) over the smoothing radius for the upper 
(lower) panels of figure §. 

The right panels of figure g show that in the presence 
of noise, the pdf's look quite different from the noise-free 
case. However the noisy pdf's still have different shapes 
from the pure noise pdf's and the cosmological models re- 
main distinguishable. The asymmetric double peak for 
the low amplitude peaks (— 2<i/<2) arises due to the 
noise maxima and minima, but it is suppressed relative to 
the pure noise case and is asymmetric due to the gravi- 
tational shear. The open and Einstein-de Sitter models 
are easily distinguishable for these low amplitude peaks, 
even though almost none of the peaks can be individually 



associated with dark matter halos. The relative number 
of very negative troughs (which are mostly noise troughs 
modulated by the fact that they are located in large voids) 
can by itself be used to discriminate models with different 
values of fl, as noted by Jain et al (f999) for the field 
pdf. As the smoothing scale is increased from 0.5' to 1', 
the signal dominates over the shape of the noise pdf, but 
sample variance becomes larger as there are fewer peaks. 
As a result, at both smoothing scales, the models can be 
distinguished at about the same level of significance. The 
error bars in the right panels are not much larger than in 
the pdf from the noise- free maps. The addition of ellip- 
ticity noise broadens the peak distribution, but the error 
bars are still dominated by sample variance in the signal. 
Similary, the primary effect of increasing a^^ is to broaden 
the pdf while not changing the error bars by much. 

It is worth noting that previous theoretical work on halo 
detection has focused on the peaks that can be individu- 
ally detected with adequate signal to noise. These would 
correspond to the parts of the pdf with i/ <^ 4 — 5, where 
sample variance is large. Clearly the bulk of the informa- 
tion on the mass function and in distinguishing models is 
at smaller or negative peak heights and can be used only 
statistically by modeling the noise pdf. 

4. DISCUSSION 

The results presented in sections 2 and 3 show that the 
peak distribution from Icnsing data has information on 
the projected mass function of dark matter halos, and is 
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Fig. 4. — Histogram of the convergence field in noise-free fields (left panel) and noisy fields (right panel). The solid lines show the open 
model and the dashed lines the EdS model. Tljup error bars are obtained from averaging over seven realizations. The field size and other 
parameters are as for the upper panels of figure 



sensitive to the cosmological model. The level of non- 
Gaussianity of the pdf is a powerful discrimant of models 
with different values of Q. Figure S shows that the mod- 
els can be distinguished from the pdf over a wide range 
of peak heights at 2 to 3-cr; by combining information at 
different peak heights and smoothing scales we can obtain 
much higher significance. Further, the third and fourth 
moments of the peak distribution for different smoothing 
scales are sensitive to the cosmological model, as expected 
qualitatively from the shapes of the distribution. We have 
also compared the peak distributions shown with a model 
with non-zero cosmological constant VIa = 0.7. The peak 
distribution in the A-model lies in-between the Einstein 
de-Sitter and open model with the same value of rJniattor- 

Beyond the dependence on the cosmological parameters, 
the peak pdf contains information on the projected mass 
function over all mass scales. It is important to test how 
accurately we can recover the pdf of peaks due to the lens- 
ing signal, and hence the projected mass function, from 
wide field lensing surveys. A straightforward approach 
is to compare the measured pdf with the predictions of 
a set of models that include the level of noise observed 
in the data. The best fit model can be found by mini- 
mizing the x^- We demonstrate in a forthcoming paper 
that the projected mass function and Q can be simultane- 
ously determined by using the normalization and shape of 
the distribution (Van Waerbeke & Jain 1999). Since we 
use information from all peak heights, not just the high-cr 
peaks that can be detected individually, the mass function 
is constrained over mass scales ranging from galactic to 
cluster sized halos. A more ambitious approach to recover 
the lensing signal would be to de-convolve the measured 
peak distribution using the analytical model for the noise. 
The nearly perfect accuracy of the analytical noise model 
(see figure H^), which we have checked for four cosmological 
models with different smoothing scales and noise distri- 
bution, gives us confidence that the lensing signal can be 
extracted from forthcoming data — either using deconvo- 
lution, or by comparing the forward convolution for a set 
of models. 

Analytical predictions of peak statistics would be valu- 
able in comparing theoretical predictions with observa- 
tions. Reblinski et al (1999) have shown that predictions 
of peak number densities based on the Press-Schechter 



model agree with the simulations for the high-u peaks. 
Detailed analytical predictions of peak number densities 
and their angular correlations by combining the Press- 
Schechter model and its extensions with our noise model 
would be useful. Further work is also needed to test the 
sensitivity of the results to the shape of the dark matter 
power spectrum. The dependence on the redshift distri- 
bution of source galaxies needs to be computed as well 
— since the level of non-Gaussianity decreases for more 
distant galaxies, increasing the redshift of source galaxies 
could mimic the effect of high-J7. 

To place the peaks approach in perspective, it is useful 
to compare it with the standard approach of measuring 
dark matter statistics using the entire field (without peak 
identification). The peak statistics rely on only a subset 
of the available information (the location and height of 
peaks, and eventually their profile), and obtaining cosmo- 
logical information from them requires additional theoret- 
ical modeling compared to field statistics. On the other 
hand, the use of peak statistics has both practical and 
theoretical advantages. Peak statistics are likely to be ro- 
bust to certain kinds of systematics — small, unknown 
errors in the galaxy ellipticities that complicate the use of 
field statistics. For example, in practice, the shear mea- 
sured from ellipticity data is multiplied by a factor larger 
than unity to account for the smearing by the point spread 
function. If this factor is estimated incorrectly, it could 
change the height but probably not the location of the 
peaks. It would then amount to a rescaling of the x-axis 
in the peak histograms, which does not change the com- 
parisons amongst the cosmological models. 

On the theoretical side, peak statistics can provide in- 
sights into the biasing of galaxies relative to the dark mat- 
ter, by allowing us to consider the two distinct components 
of biasing: first, the relation of galaxies to dark halos, and 
second, of halos to the dark matter. By combining the 
measured clustering of galaxies with that of dark halos 
measured through peak statistics from lensing data, the 
first step in the biasing of galaxies can be directly probed. 
For the second step of relating halos to the dark matter, 
we will need to use successful measurements of field lens- 
ing, or to interpret the data using theoretical models for 
the relation of halos to the dark matter. 

We have shown that the statistics of peaks provides a 



useful new approach to wide field lensing. It is complemen- 
tary to the standard statistics such as ellipticity correla- 
tions over the field, and is directly linked to the projected 
distribution of dark matter halos. The characteristic non- 
Gaussian shape of the peak distribution (its asymmetric 
double peaked shape) makes it a powerful probe of the 



cosmological model as well as a useful test of the presence 
of systematics errors. 

We are grateful to Uros Seljak, Peter Schneider, Ravi 
Sheth, Alex Szalay and Simon White for helpful discus- 
sions. We thank an anonymous referee for comments. 



REFERENCES 



Bond, J.R., Kaiser, N., Szalay, A.S., 1986, ApJ, 304, 
Narayan, R., Seitz, S., Schneider P., 1996, ApJ, 



1997, A&A 322,1 
Villumsen J.V., 1991, 



Bardeen, J.M 

15 
Bartelmann, M 

464, L115 
Bernardeau F., Van Waerbeke L., Mellier Y, 
Blandford R.D., Saust A.B., Brainerd T.G., 

MNRAS 251, 600 
Bond, J.R., Efstathiou, G., 1987, .MNRAS, 226, 
Erben, T et al., astro-ph/9907134 

Jain B., Seljak 0., 1997, ApJ 484, 560 . 

Jain, B. Seljak, U., White, S., astro-ph/9901191 

Kaiser N., 1992, ApJ 388, 272^ 

Kaiser, N. 1998, ApJ 498, 26 

Kaiser, N., Squires, G., 1993, ApJ 404, 441 



655 



Kaiser, N. et al., istro-ph/9809268 

Kruse G., Schneider R, 19y9a, MflRAS 302, 821. 

Kruse G., Schneider P., 1999b, |:istro-ph/9904192 

Miralda-Escude J. ,1991, ApJ 38U7T ^| , 

Reblinsky, K., Kruse, G., Jain, B., Schneider, P., istro-ph/9907250 

1996, MNRAS, 283, 837 

van Waerbeke L., Jain B., Kruse G. 1998, MNRAS 



1996, 



Schneider, P, 
Schneider P. 

296, 873 
Stebbins. A., 
Squires, G., Kaiser N~ 
Tyson, T., et al, 19aa 
Van Waerbeke, L., [a,stro-ph/9909160 
Van Waerbeke, L., Bernardeau, 1^'., iUellier, Y., 1999, A& A, 342, 15 
Van Waerbeke, L., Jain, B., 1999, in preparation 



astro-ph/960914g 
1996, ApJ 4V3 
in prcparatii 



65 



